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We study the thermodynamics of the universe containing perfect fluid with equation of state 
p = up, uj is variable. Here we choose uj to be a function of the red shift variable z and four different 
choices of uj has been investigated. Also the laws of thermodynamics are examined considering the 
universe bounded by the horizon (apparent) as a thermodynamical system. 

o\ '• 

O ■ I. INTRODUCTION 

o : 

In the semi-classical description of black hole physics it is found that a black hole behaves as a black 
^ body and emits thermal radiation. The temperature (known as Hawking temperature) and the entropy 
are proportional to the surface gravity at the horizon and area of the horizon [1,2] respectively. The 
\Q . Hawking temperature, entropy and mass of the black hole satisfy the first law of thermodynamics [3]. As the 
t-H temperature and entropy are determined by purely geometric quantities (namely surface gravity and horizon 
area respectively), i.e. characterized by the space-time geometry and hence by Einstein field equations, so it is 
' natural to speculate some relationships between black hole thermodynamics and Einstein equations. Jacobson 
Q H , [4] showed that Einstein equations can be derived from the first law of thermodynamics: <5Q=Tds for all local 
Rindler causal horizons with SQ and T as the energy flux and unruh temperature seen by an accelerated 
{^JQ! observer just inside the horizon. For a general static spherically symmetric space-time, Padmanabhan [5] was 
l— ~~ '■ able to derive the first law of thermodynamics on the horizon, starting from Einstein equations. 

Subsequently, this equivalence between Einstein equation and thermodynamical laws has been generalized in 
the context of cosmology. If we assume the universe as a thermodynamical system and consider at the apparent 
ON ■ horizon Ra the Hawking temperature Ta — and the entropy Sa = vR q , then it was shown that the 

first law of thermodynamics on the apparent horizon and the Friedmann equations are equivalent [6]-one can 
. ■ be derived from other. Then in higher dimensional space-time, the relation was established for gravity with 
Gauss-Bonnet term and for the Lovelock gravity theory[6,7]. As a result, it is speculated that such a deep 
relationship between the thermodynamics at the apparent horizon and the Einstein equations may give some 
clue on the properties of dark energy. 



o 
o 



X 



In the present work, we study the thermodynamics of the universe with matter in the form of perfect fluid 
having variable equation of state in form: p = uop. We consider to as a function of the red shift z in the following 
forms of two index parametrization[8]:(%) uj — ujq + uj\z (linear redshift parametrization),^,) uj = ujq + n^fj 
(Chevallier-Polarski linear parametrization) , (Hi) to = ujq + rjqr^s (Jassal-Bagla-Padmanavan parameteriza- 
tion),^^ uj = ujq + L)\z if z < 1 ,ujq + uj\ if z > 1 (Upadhya-Ishak-Steinhardt parametrization). The thermal 
quantities are expressed either as a function of volume or temperature and due to adiabatic nature of the 
thermodynamical system the entropy turns out to be a constant. Subsequently , we examine the validity of the 
thermodynamical laws for universe bounded by the horizon(event or apparent). 



II. GENERAL THERMODYNAMICAL DESCRIPTION 



Let us consider a thermodynamical system bounded in a volume V and suppose p ,p and T are the energy 
density, thermodynamical pressure and temperature of the fluid bounded by the volume. Then from the first 
law of thermodynamics [9] 



* nairwital5@gmail.com 

t schakraborty@math.jdvu.ac. in 



2 



TdS = d(pV) + P dV = d[(p + p)V] - Vdp (1) 
where S is he entropy of the thermodynamical system .Now the integrability condition [9] 

*S 8 2 S (2) 



dTdV dV&T ' 
demands the following differential relation[10] 



dp dT 



p + p T 

Combining (1) and (3) and integrating once we obtain (except for an additive constant) 

S=<^ (4) 

However, for adiabatic process entropy is constant and consequently, the first law of thermodynamics (1) 
becomes the conservation law 

d[(p + p)V] = Vdp (5) 

Note that one can obtain the relation (4) easily using the conservation relation (5)into the integrability 
condition(3). Hence for adiabatic process equation (4) may be considered as the temperature defining equation. 
If we suppose our universe to be homogeneous and isotropic FRW space-time with line element 



ds 2 = -dt 2 +a 2 {t) 



dr 2 



1 — kr 2 

then the Einsteins equations (known as Friedmann equations) and the energy conservation law are 



(6) 



2 8irGp k 

H = —o 72 ' ( ? ) 



H = 4nG(p + p) + ^ (8) 

and 

p + 3H(p + p) = 0, (9) 

where H = - is the Hubble parameter, k is the curvature scalar having values 0, ±1 for three 

dimensional space to be fiat or to have positive or negative spatial curvature and dQ,2 2 is the metric 
on unit 2-sphere.Now we consider our universe filled with perfect fluid having equation of state p = cup 
where ui is not constant but a function of the red shift z i.e. u) = w(z) .We choose the following forms 
for w(z) namely (%) uj = ujq + uj\z ,(II) uj = oj q + (f^fy, (HI) co = ojq + (IV) to = lu + u\z if 

z < 1 ,uq + u)\ if z > 1 where wo and u)\ are constants. 

Case I : ui(z) = uj + u>iz From the energy conservation equation (8) integrating once we have 

p{z) = p Q e^ lZ (l + z f 1+ ^-^ (10) 




Fig. 1(a), 1(b), 1(c) shows the variation of the energy density,temperature and velocity of sound with the red shift factor 
(given by eq.(10,ll,12)) respectively where loq > u>\ , here the values of u)q — 2 and wi = 1 . Fig. 1(d) shows the 
variation of heat capacity with red shift factor for ujq < u>i with values of uo = 1 and u>i = 3. 

Here po = 1, So = 1 



where po an integration constant. As the system is an adiabatic process so entropy is constant (So, say) and 
hence from equation (4) 



T(z) = T e 3 ^(l + z )^°-^(i +UJQ + uiz) 



with To 



j^-, a constant. The squared speed of sound v s 2 is given by 



(11) 



V S 2 {Z) 



Finally, the heat capacity is calculated to be 



dp 
dp 



u>i{l + z) 
3(1 + w) 



(12) 



C v (z) = V 



dp 35o 



dT l + z 



1 + Z l+QJ 



Limiting Values 



(i) a — > i.e. z — > oo : (from fig.l) 
p — y oo,p — > oo,T — > oo,v s — * oo and C v — » 0. 

fiij a — > oo i.e.z — > —1 : 

(a) p —t 0, p — > and T — » if ujq > cj% 

(b) p — > oo , p — > — c« and T — -> oo if w < w i < 1 + w o 



(13) 
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Fig.2(b) 



Fig.2(c) 



Fig.2(d) 



Fig. 2(a), 2(b), 2(c) and 2(d) shows the variation of the pressure, temperature ,velocity of sound and heat capacity with 
the red shift factor (given by eq.(14,15,16,17))respectively where po = 1, So = 1 and u>i < , with 

ujo = l,u>i = (--5) 



(c) p^poe 

(d) C v 



— 3wi 



So 



and p 



-p e 



if 1 + u>o = u>i and T — > Tqc 



LU - U3\ 



if wq = wi 



The limiting behaviour shows that for realistic fluid we must have ujq > u>\ . Here the energy density 
has correct behaviour, pressure is positive throughout the evolution (approaches to zero) with finite sound 
velocity. Also the third law of thermodynamics is satisfied as temperature approaches to zero with volume 
goes to infinity (see Fig. 1(b)). Further as C v > so the evolution from large temperature to T = is a 
thermodynamically stable transition without any critical point. When ujq < w\ , then pressure becomes 
negative at an intermediate stage of evolution and finally approaches to zero — a possibility of accelerating 
phase of the universe. Also as C v changes sign (shown in Fig. 1(d)) so there is a critical point indicating a 
phase transition of the thermodynamical system. For luq = u>x , the pressure cannot be negative at any stage 
of the evolution, the temperature approaches a finite non-zero value but heat capacity approaches infinite value 
as volume increases infinitely. Lastly, it should be noted that from equations (10) and (11) z can be obtained 
in principle as a function of p and T and then substituting in (10) it is possible to have p (and hence p ) 
as a function of T and hence all other thermodynamical variables can be expressed as a functions of temperature. 

Case II : w{z) = + 

Corresponding to this equation of state, p can be obtained from the energy conservation equation (9)as 



p(z) = p e 1 + z (l + z) 



3(l+wo+<*>i) 



(14) 



Due to the adiabatic nature of the thermodynamical system temperature can be obtained from equation(4) 
as a function of z namely 



T{z) = T e^ (1 + z ) 3 ("o+-!Ki (1 +UJq + z(1 +UJo+ 
The expressions for the squared speed of sound v"l and the heat capacity C v are given by 



(15) 
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and 



C v {z) = 



3S 
1 + z 



3(l + z)(l + w) 



3w 



1 + z (l + w)(l + z) 2 



Limiting Values : 



(16) 



(17) 



As a — + i.e. z — > oowc can find ( from fig. 2) p — > oo,p — > oo,T — > oo,v s 2 — ► (cj + and C v — > .In 
this limiting case oj — > o>o + w i 



Also if a — > oo i.c.z -^-lwe have w — > oo if u)\ < and w — > — oo if wi > 

— > if wi < and p — > oo if u)\ > 
(ii)p — > if wi < and p — > oo if wi > 
(Hi) T if uji < 0, T oo if wi > 
(iv)C v — ► , w s 2 — > oo if wi < and v^ 2 — > — oo if oj\ > 

wi < but w > wi then this equation of state is purely a decelerating model of the universe. But if 
ujo < |^i | then the universe changes from accelerating phase to decelerating phase. The velocity of sound has 
a finite value at early epoch and then approaches to infinite with the evolution of the universe while heat 
capacity starting from a finite value decreases to zero. The thermodynamical system obey the third law as 
T — > when volume becomes infinite. However , for oj\ > , though the universe evolve from decelerating 
phase to an accelerating phase, yet it is not a realistic model of the universe because p and T decreases in both 
the limits and square of the velocity of sound becomes negative at some intermediate instant and approaches 
— oo as volume becomes very large. 

Case-Ill: oj(z) = u>q + 



(1+z) 2 



Corresponding to this adiabatic equation of state the relevant thermodynamical quantities are given by 



p(*) = po(l + *) 3(1, ~"< 



3^(22 + 1) 2 
+^0)p 2(1 + *) 



(18) 



3uJ! (2z+l) 2 

T(z) = T (l + lo)(1 + z) 3wo e ^oW" 



(19) 



v s 2 (z) =UJQ + 



(1 + z) 2 



z + 



(1-z) 



3(1 + w) 



(20) 



C v (z) 



35 
1 + z 



wi(l - z) 



3oj 



1 + z {l+w){l + zf 



(21) 



Limiting Values : 



As a — > i.e. z — > oowe can find ( from fig. 3) p — > oo, p — > oo, T — ► oo, w s 2 — > o> an d C„ — > ^ -I 11 this 
limiting case w — ► ujq 



() 




Fig. 3(a), 3(b) shows the variation of the equation of state with the red shift factor when u\ > and 
loi < respectively. For Fig. 3(a) u?i = 2, loo = 1 and for Fig. 3(b) u>\ = (—5), coo — 1 




Fig.4(c) 



Fig.4(d) 



Fig. 4(a), 4(b), 4(c) and 4(d) shows the variation of the energy density, pressure, velocity of sound and heat capacity 
with the red shift factor (given by eq. (18, 19, 20, 21)) respectively, where po = 1, So — 1 and 
u)\ < OjWo > with loo =1,loi = (—5) 



Also if a — » oo i.c.z — > —1 we have u> — > oo if wi < and cj — » — oo if Wi > 

p — + if LUi > and p — > oo if wi < 
(ii)p — * oo if u>i < and p — > — oo if lo\ > 
fiiij T -> if wi > 0, T -> oo if wi < 
(iv)C v — > , w s 2 — > oo if cji < and u s 2 — > — oo if Wi > 

for this equation of state if u>% > then initially the universe has a decelerating phase (if too > ) and 
then it has accelerated expansion [in Fig. 3a] while for negative ui% ( with uiq > ) the universe changes from 
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decelerating phase to an accelerating phase then again there is a deceleration provided |wi| > 4|w | [in Fig. 3b]. 
Although the third law of thermodynamics is satisfied for u>i < , the thermodynamical system is not a 
realistic one as velocity of sound become imaginary at an instant. On the other hand for io\ > the density 
pressure and temperature starting from an infinite value reach a minimum value and then again increases 
unboundedly. The heat capacity starting from a finite value gradually approaches to zero with the evolution of 
the universe . Note that , although we have a deceleration— acceleration— deceleration phase for oj\ < and 
ll>o > but the only drawback is the imaginary sound speed at some finite red shift. 

case-IV : 

w(z) = w + u>iz, if z < 1 



uj +u>i if z > 1 



This equation of state represent a decelerating model of the universe throughout the evolution if loq > uj\ while 
for w < oj\ there will be a transition from deceleration to acceleration at z = — ^ . The physical parameters 
characterizing the thermodynamical system are given by 



p(z)^ Po e^ z (l + zf [1+UJa - u)l) ifz<\ and p (l + z) 



3(w +wi) 



if Z> 1 



(22) 



T(z) = T e 3 ^(l + z) 



(oio-oii) 



(l+w +wiz) ifz<\ and T (l + z) 



3(w +^i) 



if z > 1 (23) 



9 , s (1 + Z) 



aw 

Limiting Values : 



35o 
1 + 2 



3(1+ a;) 
3w 



and 



1 + 2! 1 + W 



z/z < 1 



and 



w + Wi if z > 1 
So 



l + z d 



i/z>l 



(24) 
(25) 



For the volume of the universe approaches zero or infinity the limiting value of these thermodynamical 
parameters arc given by as follows: 

(I) , a — > i.e. z — > oo: 

— > oo, p — > oo, T — > oo if (w + > 
(ii)p 0, p -> 0, T if (w + wi) < 
(raj/o -> po, P -> 0, T ->■ po if (w + = 
(iv)v s 2 — > (wo + w i) ! Cd — > and w — > wo + wi 

(II) . a — > oo i.c.z — > —1 : 
(%) p, p, T — > if uji < uj 

(ii) p, p, T — > oo if w < Wi < 1 + wo 

(iii) v 3 2 



{ujo-uji) , C„ -» and 



w — > Wq — 



It is clear that wo < Wi is not physically reasonable as velocity of sound become imaginary for infinitely 
large volume. For ojo > oji the third law of thermodynamics is satisfied, velocity of sound has finite non zero 
value and heat capacity gradually increases to a finite value throughout the evolution. 
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III. VALIDITY OF LAWS OF THERMODYNAMICS 

Suppose the line element (see equation(6)) for FRW model is written in the form 

ds 2 = h ab dx a dx b + P 2 dn 2 2 (26) 

where f = ar is the area radius (or geometrical radius) and h ab — diag(—l, ^^ r ^ ) is the metric on the 

2 dimensional hyper surface (x° = t, x 1 = r ) .Then the dynamical apparent horizon which is a marginally 
trapped surface with vanishing expansion, is given by the relation 

h ab d a fd b f = (27) 

i.e. 

/ k 2 

.2 / rr2 



1 - Ra'^H" + - j = 

The above relation shows that for flat space(fc = 0) -fi^the radius of the apparent horizon coincide with the 
Hubble horizon. The cosmological event horizon is defined as 

f°° dt [ x da 

R ° = "l T-J. BP (28) 

The thermodynamical quantities namely entropy and temperature have the following expressions at the 
horizon 

A 7ri?/ 2 1 , . 

s < = w = -(r = (29) 

where / = A or E according as we consider apparent or cosmological event horizon. 

The first law of thermodynamics on a horizon (apparent or cosmological event horizon) can be written as 

- dEj = TjdSi (30) 

Where —dEj is the amount of energy crossing the horizon during the time interval dt .For the present 4D 
FRW space time the expression for dEi is given by 

- dEi = ATrR I 3 H(p + p)dt , (31) 

where p and p are the energy density and thermodynamic pressure of the matter in the universe bounded by 
the horizon and the equation of state is given by (as in the previous section ) p=uop where u is a function of 
the red shift variable z. Now for the horizon Ri the change of entropy is given by (choosing G=l) 

dS7 = 2irR I dR I (32) 

For the apparent horizon Ra using the definition and the Friedmann equation (7) we have from above (on 
simplification) 

dS A 2 = 8tt 2 R A i H(p + p)dt (33) 

and hence 

T A dS A = 4ttR a 3 H( P + p)dt 
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So from equation (31) we have 

-dE A = T A dS A 

Thus the first law of thermodynamics is obeyed at the apparent horizon. The validity of this law does not 
depend on the equation of state of the fluid in the universe bounded by the apparent horizon. 
On the other hand, the radius of the event horizon cannot be obtained in a close form (only in the integral 
form given by equation (28) for arbitrary equation of state and as a result the right hand side of equation (7) 
can no longer be simplified further . Hence no conclusion can be made about the validity of the first law of 
thermodynamics at the cosmological event horizon. 

To examine the validity of the second law at the horizon (apparent or event) we first study the entropy enveloped 
by the horizon. Using the Gibb's equation [11] one can relate the entropy of the universe inside the horizon to 
its energy and pressure in the horizon as 

TdS m = dEm + pdV (34) 

Where V = is the volume bounded by the horizon, E m = 47rJ ^ J p is the energy inside the horizon and 

S m is the entropy of the matter distribution inside the horizon. We now consider the apparent horizon as 
the boundary of the universe. For thermodynamical equilibrium we can choose the temperature T A on the 
horizon as the temperatue of the matter inside R A . Then from the above Gibb's law(i.e. equation (34)) [12] 

^=2^(1 + 9) (35) 

In deriving this equation we have used the Friedmann equation (8) and energy conservation relation(9).Hcrc 
q = —1 — jp- is the usual deceleration parameter. In the similar way from equation (33) we have 

^ = 2*R A (l + q) (36) 

and hence 

d{Sm dt = 27tRa{1 + q)2 ~ ° (37) 

Thus the total entropy of the matter distribution inside the (apparent) horizon and the entropy of the 
horizon always increases with the evolution of the universe i.e. the generalized second law of thermodynamics 
is valid for the Universe bounded by the apparent horizon. However, the entropy of the matter distribution 
inside the apparent horizon increases for decelerating phase of the universe while the entropy decreases in the 
accelerating phase. 

As we have mentioned earlier that there is no explicit analytic form of Re so it is not possible to obtain an 
analytic expression for both ^J^- and ^f- and hence we cannot draw any conclusion about increase or 
decrease of the entropy with the evolution of the universe. So no positive conclusion is possible for the validity 
of the thermodynamical laws at the cosmological event horizon. 



IV. DISCUSSIONS 

In the present paper we consider the general thermodynamics of the universe filled with perfect fluid having 
equation of state p = cop. Here uj is chosen as a function of the red shift variable z . Four different choices for 
lo are taken and thermodynamics with asymptotic limits have been discussed. These choices of u> are recently 
shown to be in good agreement with current observations in different ranges of z .The third choice of w namely 
that of Jassal-Bagala-Padmanavan seems to be interesting because it shows a transition from deceleration to 
acceleration and then again deceleration as expected from observation. 

Next section deals with validity of the thermodynamical laws for the universe filled with perfect fluid ( p = u>p, u> 
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is variable) and is bounded by apparent or cosmological event horizon. As the event horizon radius can not 
be evaluated in closed form (only in an integral form) so it is not possible to infer about the validity of the 
thermodynamical laws. However for the apparent horizon, both first and second law of thermodynamics are 
obeyed for general oj (without specifying it). For the validity of the second law we have started with the Gibb's 
equation and calculated the rate of change of the sum total of the entropy of the matter and the entropy of the 
boundary which is always positive. This is known as generalized second law of thermodynamics. For future work 
it will be interesting to consider holographic model of dark energy so that an analytic expression for the radius 
of the event horizon can be obtained and then it will be possible to examine the validity of thermodynamical 
laws at the event horizon. 
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